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Abstract 

We discuss the role of lightest neutrino mass (mo) in the neutrino mass matrix, defined in a flavor basis, 
through a bottom-up approach using the current neutrino oscillation data. We find that if mo < 10 _3 eV, 
then the deviation 5M U in the neutrino mass matrix from a tree-level, say tribimaximal neutrino mass 
matrix, does not depend on niQ. As a result 5M„s are exactly predicted in terms of the experimentally 
determined quantities such as solar and atmospheric mass squared differences and the mixing angles. On 
the other hand for mo ^ 10~ 3 eV, 5M U strongly depends on mo and hence can not be determined within 
the knowledge of oscillation parameters alone. In this limit, we provide an exponential parameterization 
for bM v for all values of mo such that it can factorize the mo dependency of 5M V from rest of the oscillation 
parameters. This helps us in finding 5M V as a function of the solar and atmospheric mass squared differences 
and the mixing angles for any values of mo . We use this information to build up a model of neutrino masses 
and mixings in a top-down scenario which can predict large #13 perturbatively. 
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I. INTRODUCTION 



Compelling evidences of neutrino oscillations observed in solar, atmospheric, and reactor exper- 
iments indicate that neutrinos are massive and hence they mix with each other In the basis 
where the charged lepton masses are real and diagonal, the mixing matrix is given by 



PMNS 



iScp \ 



Cl2 C13 Si 2 C13 S13 C 

— C23 S12 — S23 C12 S13 e lS ° p C23 C12 — S23 S12 S13 e l6cp S23 C13 

\ ^23 s 12 — C23 C12 S13 e tScp — S23 C12 — C23 S12 S13 e lScp C23 C13 J 



where c 



cos 6^, Sjj = sin^jj and #cp is the Dirac CP violating phase which appear in oscillation 



experiments. The diagonal matrix P = diag(e J9i:L , e"^ 2 , 1), consists of two Majorana phases <f>\ and 
(p2 which are not relevant in neutrino oscillation experiments (3), Q|. However, they affect lepton 
number violating amplitudes such as neutrinoless double beta decay. Thus in a flavor basis, where 
charged leptons are real and diagonal, the neutrino mass matrix is given by: 



M v = Up MNS diag( 



(2) 



where mi, m 2 , m 3 are the mass eigenvalues. 

In the last decade, data from solar and reactor neutrino experiments have provided information 
on the sign and magnitude of Amj| and a precise value of 9u |s, 6|. The atmospheric parameters 



I Am^ tm | and 623 have been measured and their precision will be increased by T2K [7] and NOz/A [8]. 
The sign of solar mass splitting Ar/ig is precisely known, while the sign of atmospheric mass 
splitting Am^ tm is still unknown. This opens up a possibility of whether neutrino masses follow 
normal ordering, i.e., mi < 7712 < or inverted ordering, i.e., < m\ < 1712. In other words, 
the lightest mass, either mi (normal ordering) or (inverted ordering) is yet to be determined. 



Recent measurement from T2K 



MINOS 



lfj, Double Chooz Daya Bay [12j, and RENO [13 1 
confirms a non-zero value of #13 at 5a confidence level. This opens up a range of possibilities to 
measure the sign of the atmospheric mass splitting and the unknown CP phase Sep- 

The absolute mass scale of neutrino is hitherto not known and can only be measured in a tritium 
beta decay experiment. The KATRIN experiment, which will investigate the kinematics of tritium 
beta decay, aims to measure the neutrino mass with a sensitivity of 0.2 eV [ijj]. At present the 
best upper limit at 95% confidence level on the sum of the neutrino masses comes from the cosmic 
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microwave background data and is given by 151 ] : 



Y, mi < OMeV . (3) 

i 

Once the absolute mass scale of the lightest neutrino, either mi in normal ordering or m 3 in inverted 
ordering, is determined one can reconstruct the neutrino mass matrix in the flavor basis using the 
experimental values of the elements of £/pmns matrix. This may unravel exact symmetries in the 
neutrino mass matrix. 

In this paper we develop a perturbative bottom-up approach to unravel the flavor structure of 
neutrino mass matrix. We set the full neutrino mass matrix: M v = (M„)o + 5M U , where 8M V is 
the perturbation around the tree-level mass matrix (M u ) , which is determined using some of the 



well known mixing scenarios such as tribimaximal (TBM) mixing [if]], bimaximal (BM) mixing 17] 
and/or democratic (DC) mixing isj]. Among all these, the TBM is closer to the experimentally 
observed mixing pattern and hence mostly studied. All these mixing scenarios, however, predict 
$13 to be zero and hence ruled out by the recent measurement on the reactor neutrino angle. 
However, a perturbative approach to realize a large value of #13 is still a viable option. So, we 
assume that the non-zero value of #13 is generated perturbatively. Using the 3a range of values 
of the elements of C/pmns matrix we determine 8M V as a function of the lightest neutrino mass, 
say m , where m = ra\ in the normal ordering and m = m 3 in the inverted ordering. In this 
way we find that for mo < 10 -3 eV, all the 5M V are independent of mo and hence the lightest 
neutrino mass does not play any role in the perturbative determination of neutrino mass matrix. 
Such models, for example, can be generated in two right-handed neutrino extensions of the standard 
model (SM). However, for mo £ 10~ 3 eV the lightest neutrino mass plays an important role in the 
perturbative determination of neutrino mass matrix. We factorize the mo dependency of 5M V using 
an exponential parameterization: 5M v (mo) oc <5M„| mo= oExp[— mo/0.1 eV]. This helps us in finding 
the required perturbations of 5M U as a function of experimentally determined quantities such as 
solar and atmospheric mass squared differences and the mixing angles. Using this results from 
the bottom-up approach we provide an example in the top-down scenario which predict large 6*13 
perturbatively. 

The paper is organized as follows. In section Hi Al we start with a brief description of the bottom- 
up approach that is developed in this paper to find the desired perturbation of the mass matrix in 
the flavor basis. Then, in section III B\ we use the knowledge of our current experimental results to 
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find the deviation in the mass matrix in the flavor basis using the tree level mass matrix that is 
generated using the TBM mixing ansatz. In section UTTl we suggest a typical top-down approach in 
which the perturbation to the TBM mixing is obtained. We then fit the model parameters using 
the data from the bottom-up approach of section III Bl and conclude in section IIVI 



II. PHENOMENOLOGY 



A. Bottom up approach 



The recent global fit to the neutrino oscillation data has ruled out the possibility of a zero reactor 



angle at 10.2a confidence level 19, 20 1 . Evidence of non zero #13 at 3a level was first established 

n nr n n 

by data from T2K [7], MINOS [10| and Double Chooz [ljj. More recently, Daya Bay [12| and 
RENO |L3J] experiments confirm large #13 at more than 5a confidence level (CL) from the reactor 
v e — > v e oscillations. The current best fit values and the 3a allowed values of all the mixing 
parameters are summarized in Table. [B The values written within brackets are for inverted ordering 



Oscillation parameters 


Best fit value 


3a range 


Am| x (10 5 eV 2 ) 
|Am 2 tm |x(10 3 eV 2 ) 
sin 2 6>i2 
sin 2 6*23 
sin 2 6>i3 


7.62 
2.55 (2.43) 
0.320 
0.613 (0.600) 
0.0246 (0.0250) 


[7.12 - 8.20] 
[2.31 - 2.74] ([2.21 - 2.64]) 

[0.27 - 0.37] 
[0.36 - 0.68] ([0.37- 0.67]) 
[0.017 - 0.033] 



TABLE I: Current status of oscillation parameters taken from Ref. 
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of the neutrino mass spectrum. At this juncture we note that the only unknown quantity in the 
neutrino mass matrix is the lightest neutrino mass apart from the CP violating phases (one Dirac 
and two Majorana phases) 

The neutrino mixing matrix derived from the current best fit parameters is 



U EXPT 



0.81908 0.54773 0.17058 



-0.49795 0.53116 0.68551 
^ 0.28487 -0.64642 0.70780 J 



(4) 
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and the corresponding 3a allowed values of the mixing matrix is 

' 0.78 - 0.85 0.51 -0.60 0.13 - 0.18 N 
%xpt = -(0.39 - 0.57) 0.36 - 0.64 0.59 - 0.82 , (5) 
v 0.19 - 0.44 -(0.54 - 0.76) 0.56 - 0.79 , 

where we have assumed 5c p to be zero for simplicity. The corresponding neutrino mass matrix in 
the flavor basis can be written as 

(M^expt = Uexpt diag(mi, m 2 , m 3 ) f/| XPT (6) 

where (M^)expt is a real symmetric matrix. Therefore, it is described by three masses: mi, m2, 
and three mixing angles: #12, #23, #13- The mass eigenvalues mi, m2, in the basis where 



the charged leptons are real and diagonal can be expressed as: m2 = J ArriQ + mf and m 3 



Am^tm + m i f° r normal hierarchy and mi = y mj + Am^m and m2 = J m| + Am^m + Am© for 
inverted hierarchy. 

We define the perturbed neutrino mass matrix 5M U as 

5M U = (M,) EXP T - {M v ) , (7) 

where (Mj,)expt is the neutrino mass matrix given by Eq. [6] and (M u )q is the neutrino mass matrix 
corresponding to a mixing matrix Uq, that is obtained in a specific model such as TBM mixing, 
BM mixing and/or DC mixing etc. This 5M U is a symmetric matrix whose elements denote the 
amount of perturbation that is needed to add to the (M V ) Q so that it is consistent within the 3a of 
the experimental results. In order to gauge the size of the SM U matrix, we do a random scan of all 
the parameters: #13, #23, 612, Amg and |Am^ tm | within their 3a allowed range of values. Note that 
(Mj,) depends only on the value of the lightest neutrino mass m as the mixing angles are expected 
to be fixed via certain symmetries as done in TBM, BM and/or DC mixing scenarios. Hence, to 
see the effect of m on the elements of 5M U we vary m in the range [10~ 9 , 10] eV. We then define 
the predicted neutrino mass matrix (M„)m of any model as: 

(M U ) M = (M U ) + 5M U . (8) 

Note that the matrix elements of (M u )m and 5M U are function of the oscillation parameters #12, #23, 
and #13 as well as of the lightest neutrino mass mo- We then perform a naive statistical analysis. 
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We compare the elements of the (M u )m matrix to the experimental data of Eq. [6] by a x 2 function 
which is defined by: 



X 



{Ml 



IEXPT 



-{M V ) M /m a 2 = 5M C V -5M V /m a 



(9) 



where a = J (M£) EXPT and 5M^ = (M^) EX pt — (^C)o- Note that (M^)expt is the experimental 
mass matrix corresponding to the best fit values of the oscillation parameters and is diagonalized 
by the mixing matrix of Eq. HI It is worth noting that as SM P — > 5M£, we get the minimum \ 2 and 
hence 6M£ is the required amount of perturbation that should be added to (M v )q. We will discuss 
the properties of the 5M U matrix using the TBM mixing ansatz. Although we focus mainly on the 
TBM mixing, some of the features presented here are applicable to other tree-level mass matrix. 



B. Perturbation in Tribimaximal mixing Scenario 



We proceed to discuss the tribimaximal (TBM) mixing scenario. For the TBM mixing , we have 
S12 = l/v3, S23 = l/\/2, and S13 = 0. We know that apart from the recently measured value of 
#13, the solar and atmospheric mixing angles are in a very good agreement with the TBM ansatz. 
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Modifications to this TBM scenario have been studied by many authors 

We begin by writing the TBM mixing matrix in the standard parameterization of Eq. dj that is 

/ . /o/q . FTR n \ / nsifi5n 057735 n \ 



TBM 



2/3 ^1/3 
^176 /l73 ^172 



0.81650 0.57735 
-0.40825 0.57735 0.70711 
^ 0.40825 -0.57735 0.70711 J 



(10) 



The corresponding neutrino mass matrix in the flavor basis and the perturbed matrix 5M U are 



(M u )tbm = C/TBMdiag(mi, m 2 , m 3 )U^ BM 



6M U = (M v ] 



EXPT 



TBM- 



We wish to determine the elements of this SM U matrix in a bottom-up approach. We follow the 
procedure given in section Hi Al and perform a random scan over all the oscillation parameters within 
their 3a range of values. In Fig. (TJ we show 8M v {i,j) versus \ 2 f° r m o — 10~ 9 eV. The minimum 
of each curve corresponds to the value of 5M°(i,j). A similar plot is shown in Fig. [2] for inverted 
hierarchy as well. We know that the neutrino mass matrix in the flavor basis does depend on 
the value of the lightest eigenvalue m . In order to see the effect of m on the value of 5M£(i,j) 



qualitatively, we vary m in the range (10 , 10) eV while keeping other parameters such as #12, 



13, 




FIG. 1: x 2 °f SM u (i,j) for normal ordering of neutrino masses with TBM mixing in tree-level. 

623, \Aml tm \, and Amg fixed at their best fit values. In Fig. [31 we plot 5M£(i,j) versus mo and it is 
clear that for mo < 10~ 3 eV, the value of 8M°(i,j) remains constant. In other words, SM u (i,j) are 
independent of the value of m . However, for m > lCT 3 eV, the value of SM u (i,j) keeps changing. 
In the degenerate limit when mi ~ m2 ~ m.3, 5M u (i,j) become zero. This particular feature is 
common to both normal and inverted hierarchies as can be seen from Fig. [3j To appreciate the 
plots in Fig. [HI we write down 5M u (i,j) as a function of neutrino masses and mixing angles: 

5M U {1, 1) = ((c 12 c 13 ) 2 - 2/3)771! + ((s 12 c 13 ) 2 - l/3)m 2 + s 2 l3 m 3 
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FIG. 2: x 2 °f 8M„(i,j) for inverted ordering of neutrino masses with TBM mixing in tree-level. 

5M V (1, 2) = (ci2Ci 3 (-si 2 C23 - C12S23S13) + l/3)mi 

+ 0l2Cl3(Ci2C 23 - S12S23S13) - l/3)m 2 + Si 3 Ci3S2 3 m3 

<JAf„(l,3) = (ci 2 ci 3 (si2S23 - C12C23S13) - l/3)mi 

+ (Sl2Ci 3 (-Ci2S 2 3 - S12C23S13) + l/3)m 2 + Si 3 Ci3C23m 3 

SM V (2,2) = ((-S12C23 - C12S23S13) 2 - l/6)mi 

+ ((C12C23 - S12S23S13) 2 - l/3)m 2 + ((C13S23) 2 - l/2)m 3 
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1 o" 4 1 cr 3 1 o"' 

m„(eV) 




FIG. 3: Variation of SM V with respect to the lightest neutrino mass for NH (left) and IH (right) of neutrino 
masses with TBM mixing. 

8M U (2, 3) = ((-s 12 c 23 - C12S23S13X.S12.S23 - C12C23S13) + l/6)mi 

+ ((C12C23 - si 2 S23Si 3 )(-ci2S23 - S12C23.S13) + l/3)m 2 + (S23C13C23C13 - l/2)m 3 

<yM„(3,3) = ((S12S23 - C12C23S13) 2 - 1/6)77?,! 

+ ((-C12S23 - S12C23S13) 2 - l/3)m 2 + ((ci 3 c 23 ) 2 - l/2)m 3 . (12) 

Note that the higher eigenvalues 7772 and 7773(7771) can be re-expressed as a function of lightest 
neutrino mass mo = 7771(7773). Therefore, it is obvious that all the elements: 8M v (i,j) are function 
of the only unknown quantity m . From Fig.[3]we see that 5M V (1, 1), SM„(1, 2), 5M V (1, 3), 6M U (2, 2) 
start with a positive value, while SM U (2,3) and SM U (3, 3) start with a negative value. An exactly 
opposite spectrum is observed in case of inverted hierarchy as expected. In either case we observe 
that for mo < 10~ 3 eV, SM u (i,j) are independent of mo- Therefore, in this limit it is reasonable to 



set mo ~ 0. In the opposite limit when mo >> y|Amat m | = 4.9 x 10~ 2 eV, we get 7773 « 7772 ~ mi. 
In this limit M v (i,j) w mo (/), where I is the identity matrix. Hence all the 5M u (i,j) are zero. For 
m ~ 10~ 3 eV, it is comparable to solar and atmospheric values. Therefore from Fig. [3] we see an 
appreciable effect of m on various 5M u (i,j). We have factored out this dependency of all 8M u (i,j) 
on mo using an exponential parameterization: 



6M„(m , yjAml, \/\Aml tm \, 9 12 , 9 13 , 623) = 5M u (m = 0, \f Am%, \/\Aml tm \, 9 12 , 613, 623) 

xExp(—m /a). (13) 

where a is determined from a fit to the recent experimental data and found to be a ~ 0.1 eV. As a 
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result we could do all our analysis for mo < 10 _3 eV which is equivalent to setting mo = 0. Then 
we generalize it to any value of mo using Eq. [T3j 

We first compare the perturbed matrix elements 8M£(i,j) with the elements of tree level mass 
matrix (M i/ )xbm(^ j)- in case of normal hierarchy, the experimental mass matrix, the tree-level 
neutrino mass matrix, and its deviations from the central values are given by: 



{Ml 



lEXPT 



10~ 2 eV 



0.397 0.829 0.145 



0.829 3.191 2.128 
^ 0.145 2.128 2.335 J 



( 



TBM 



10- 2 eV 



0.291 0.291 -0.291 



0.291 2.816 2.234 
y -0.291 2.234 2.816 J 



(14) 



and 



8Mf. 



10- 2 eV 



( 0.106 0.538 0.436 N 

0.538 0.375 -0.106 
\ 0.436 -0.106 -0.481 J 



(15) 



Similarly, for inverted hierarchy, the mass matrix (M^)expt and (Mj,)tbm are, 



EXPT 



10-2 eV 



( 4.830 -0.577 
-0.577 2.061 

V 



0.517 \ 
2.379 

0.517 -2.379 3.044 



(M, 



t/jTBM 



10-2 eV 



4.956 0.026 -0.026 
0.026 2.490 -2.490 
-0.026 -2.490 2.490 



(16) 



and the corresponding perturbed matrix is 



SMI 
10-2 eV 



We introduce a new parameter e as: 



I -0.125 -0.602 -0.491 N 
-0.602 -0.429 0.112 
-0.491 0.112 0.554 



(17) 



18) 



(M„) TBU (i,j) 

such that it will give a measure of required perturbation with respect to the corresponding tree-level 
value. Thus for normal and inverted hierarchies of neutrino mass spectrum, we get: 



(e c ) NH 



I 0.364 1.848 -1.848"^ 

1.848 0.133 -0.047 
^ -1.848 -0.047 0.171 J 



IH 



( 0.02 23.154 18.885 N 
23.154 0.172 0.045 
18.885 0.045 0.222 



(19) 
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From Eq. (|T9|) . it is clear that (1,2), and (1,3) elements of the TBM mass matrix needs to be 
modified largely to be consistent with the experiment. The perturbations of rest of the elements 
are small in comparison to their tree-level masses. In case of NH the modification is mild, while it 
is significant in case of IH case. 

We now try to explore the #13, #23, and 612 dependency of all the elements of the neutrino mass 
matrix in the flavor basis where the charged lepton is assumed to be diagonal. In order to see 
the effect of a particular oscillation parameter on the value of SM U , we vary that parameter in the 
3cr allowed range and do a marginalization over the other oscillation parameters. From Eq. (Tl2l) 
we see that in either case of NH and IH the dependency of 5M U (1,2) and 5^^,(1,3) on 9 12 is 
almost negligible because of the cancellation in the second term of each expression. However, both 
8M V (1, 2) and 5M U (1, 3) do depend on 2 3 and #13 in either case of NH and IH as can be seen from 
Fig. HI The dependency of 6M U (1, 2) and <5M„(1, 3) elements on #13 and 623 can be understood from 
the following analytical approximation. In case of NH, m x — > 0. Hence we have m 3 = \J Am\ tia 
and m 2 = yj Am®. As a result the second term of 5M U (1,2) and 5M u (l,3) is proportional to 
Attiq <C \J Am 2 tm . Moreover, there is a cancellation in the second term for which the second term 
is suppressed in comparison to the third term. Therefore, we get 

<JM„(1,2) w \j Am 2 atm sin 9 13 cos 6>i 3 sin 9 23 

<JM„(1,3) w J Aml tm sin 9 13 cos 9 13 cos 6 23 . (20) 
On the other hand in case of IH, 777.3 — ► 0. Hence we have 7772 ~ mi = J Am^. As a result we get 

|5M„(1,2)| \J Aml tm sin 9 13 cos 9 13 sin 6» 23 

1^(1,3)1 « y/ Aml tm sin 9 13 cos 9 13 cos 9 23 . (21) 

Thus in either cases the dependency of 8M U (1,2) and 5M V (1,3) elements on 13 and 9 2 3 are very 
similar. This can be checked from Fig. HJ It is evident from Fig. Hlthat 5M U (2,2), 8M U (2,3), and 
5M„(3, 3) elements does not depend on # 13 and 9 12 at all. It only depends on the value of #23- To 
understand this feature quantitatively, we write down the matrix elements SM U (2,2), <5M„(2,3), 
and 5M U (3, 3) for m = 0. For normal ordering of the neutrino mass spectrum, we have 

6M V (2,2) « (cos 2 fl 13 sin 2 9 23 - \) ^ Aml tm , 

5M U {2,3) pa (sin 023 cos0 23 cos 2 13 - -) \/Aml tia , 

5M U (3, 3) pa (cos 2 9 13 cos 2 23 - \) ^ Amt tm (22) 
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0.018 0.02 0.022 0.024 0.026 0.028 0.03 0.032 
sin 2 e 1 o 



0.018 0.02 0.022 0.024 0.026 0.028 0.03 
sin 2 e 1 o 





0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 0.35 0.36 0.37 
sin 2 e 1P 



0.008 
0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 




(1,1) 




(1,2) - 


IH ■ 


' (1.3) 


(2,2) 




■ (2,3) 




(3,3) 







0.27 0.28 0.29 0.3 0.31 0.32 0.33 0.34 0.35 
sin 2 e 1P 



FIG. 4: 5M u (i,j) as function of mixing angles for tuq = in case of NH (left) and IH (right) of neutrino 
mass matrix. 

Clearly, SM U (2, 2), 5M V (2, 3), and 6M U (3, 3) does not depend on 12 at all. There is #i 3 dependency 
and is of similar order if 823 is near to the TBM value. However, once the value of #23 deviates 
away from the 823 (TBM) value, all these matrix elements mainly depend on the value of 823 and 
the value of 5M u (i,j) increases as we move away from the TBM value. In case of IH, 7713 — > and 
a similar pattern for all the SM u (i,j) is expected. This can be easily read from the right pannel of 
Fig.H 
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III. TOP-DOWN APPROACH 



Now that we know the perturbed matrix from the bottom-up approach, we wish to construct 
models of neutrino masses and mixings what we call top-down approach. Our main objective here 
is that we will use the predictions of bottom-up approach as guide lines for building neutrino mass 
matrix whose tree-level mixing is governed by a TBM ansatz. We modify the TBM mixing matrix 
using some perturbation so that the resulting values are consistent with the results that are obtained 
using the bottom up approach of section. Ill B[ We begin by proposing a model where the neutrino 
mixing is described by Um = ^tbm V, where V is the perturbed mixing matrix around the TBM 
value and in general is given by: 



V = V(l,2)l/(l,3)l/(2,3), 



(23) 



V = V(1,2)V(1,3) V(2,3) 



with V(i, j) being given by an orthogonal rotation matrix in the plane of neutrino mass matrix. 
Let us write down V(i,j) explicitly such that the perturbed mixing matrix is 

I 1 P o\ ( 1 a\ ( 1 

-0 10 10 1 7 
\ 1 J \ -a 1 J \0 -j 1 J 

' 1 /3 — cry a + /?7 ^ 
= -P 1 1-pa , (24) 
\-a -7 1 J 

where the determinant of V(i,j) matrix is assumed to be unity. Since Um is the diagonalising 
matrix of the proposed neutrino mass matrix, we get 

= U T M (M U ) M U m = U T M ((M,)tbm + 5M V ) U M 

= V T ^tbm ((M„) T bm + 5M V ) U TBM V , (25) 

where, (M u )m is the mass matrix in a flavor basis where charged leptons are real and diagonal. We 
now define a matrix X such that 



X = VM^V 1 =U* BU ((M V ) TBM + 5M„)U TBU , 
which leads to the following equations in terms of a, P, and 7 in the mi — > limit. 

In = m 2 p 2 + m 3 a 2 , 



(26) 
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X 12 = m 2 (3 + (m 3 - m 2 )a 7 , 
X 1Z = m 3 a +{m 3 - m 2 )/?7 , 
X 22 = m 2 + m 3 7 2 , 
X23 = -m 3 a /3 + (m 3 - m 2 )7 , 

^33 = "7.3 + m 2 7 2 , (27) 



where m 2 = yArriQ and m 3 = yAm 2 tm , respectively. Here we have ignored higher order terms in 
a, (3, and 7. In Table. [TT1 we report all the elements of the matrix X obtained using the bottom up 
approach for the TBM mixing scenarios. We need to find the range of a, (3, and 7 so that X(i,j) 



TBM 


X(l,l) 


X(l,2) 


X(l,3) 


X(2,2) 


X(2,3) 


X(3,3) 


Central values 
3cr range 


0.0203 
[0.87 x 10~ 4 , 0.269] 


0.0487 
[-0.0937,0.118] 


0.315 
[0.0141,1.15] 


1.01 
[0.839,1.17] 


0.747 
[-0.153,1.05] 


4.89 
[4.57,5.14] 



TABLE II: Values of X(i,j) x 10 2 for TBM mixing from the bottom up approach. 



values computed from Eq. [27] are consistent with the values reported in Table. [TTJ We get the range 
of a, (3, and 7 by fitting AT(1,2), ATI, 3) and X(2,2) in Eq. [27] with the corresponding ATI, 2), 
X(l, 3), and X(2, 2) ranges from Table. HI1 Once the range in a, (3, and 7 are known, we then plug 
those in Eq. [27] to get the range in X(l, 1), A"(2, 3), and X(3, 3) and are shown in Fig. [5j In order 
to realize the dependence of Eq. [27] on number of parameters, we set 7 equal to zero and compare 
it with 7 7^ case. The results are shown in Fig. [6] It is clear that once we set 7 = 0, we do not 
get a very good fit for X(2, 2) and X(2, 3). We did the same exercise by setting (3 = 0. We fix the 
range of a and 7 from X(l, 3) and X(2, 2). Once the range of a and 7 are fixed, we get the range 
of X(l, 1), X(l, 2), X(2, 3), and X(3, 3) using Eq. [271 and is shown in Fig. [7J Similarly, once we set 
a = and fix the range of (3 and 7 from X(l, 2) and X(2, 2), we get the range of X(l, 1), X(l, 3), 
X(2, 3), and X(3, 3) using Eq. [27J The resulting ranges are shown in Fig. [8] It is clear from Fig. [5] 
that, we get a very good fit with the experimental data set if all the 3 model parameters (a, (3, 7) 
are non zero. However, setting f3 and 7 to zero also gives a quite good fit to the experimental 
data even though there are some discrepancies between the data and the theory for some matrix 
elements. But it is evident from Fig. [S]that, non zero values of j3 and 7 alone can not reproduce the 
experimental results. Thus, to have a reasonably good fit we at least need to have a non zero a, i.e.. 
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FIG. 5: Range of X(i,j) as a function of a for a, f3, 7 7^ 0. The theory ranges are denoted by red whereas 
the experimental ranges from the bottom up approach are denoted by green. 

the perturbation in the (1, 3) plane. Again, since we can set both (3 and 7 to be zero, the order in 
which we do the rotations to construct the perturbed matrix V is not very important. We further 
note that if V is product of atleast two matrices (that is V is parameterised by atleast two variables) 
then the parameters in ^-matrix can not be considered as true-perturbation around their tree-level 
values. Therefore, the new parameters a,/3,j... are arbitrary. However, if V is parameterised by a 
single variable, say a, then it can be considered as a true-perturbation around its tree-level value. 
Such a case can be obtained using an A4 symmetry in a type-II seesaw framework 16|. 



IV. CONCLUSION 

The large value of #13 predicted by Daya Bay and RENO put a stringent constraint on theoretical 
model buildings for neutrino mixing. All the well known mixing ansatz such as TBM mixing, BM 
mixing, and DC mixing predicts #13 to be zero and hence not consistent with experiment. However, 
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FIG. 6: Range of X(i, j) as a function of a with 7 = 0. The theory ranges are denoted by red whereas the 
experimental ranges from the bottom up approach are denoted by green. 

a lot of phenomenological models have been proposed in order to explain the non zero #13 where 
a large non zero value of the third mixing angle #13 is generated through a perturbation to these 
mixing scenarios. In this context, we use a perturbative bottom up approach where the neutrino 
mass matrix in the flavor basis derived from the current experimental 3a allowed ranges is written 
as a perturbation around a tree level mass matrix that is determined using the TBM mixing ansatz. 
However, our analysis is more general and can be applied to any other ansatzs which predicts #13 = 
at the tree-level. 

We know that neutrino mass matrix in the flavor basis does depend on the lightest neutrino 
mass (mo) in addition to the other oscillation parameters. In this context it is worth mentioning 
that for m <C 10~ 3 eV, the perturbed elements does not depend on the value of m . However, 
in the opposite limit, where m ^> 10 _3 eV, we find that the perturbed elements have exponential 
dependency on the value of mo- We factorise the mo dependency of all 5M u (i, j) using an exponential 
parameterization as given in Eq. ( |T3|) . This helped us in determining the exact dependency of the 
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FIG. 7: Range of X(i,j) as a function of a with j3 = 0. The theory ranges are denoted by red whereas the 
experimental ranges from the bottom up approach are denoted by green. 

perturbed matrix elements on the value of #13, #23 and #12. In order to gauge the size of the 
perturbation to each element of the mass matrix, we first compared the perturbed matrix elements 
with the corresponding tree level mass matrix derived from the TBM mixing ansatz and find that 
only the (1, 2) and (1, 3) elements needs to be modified to be consistent with the experimental data. 
We also show the exact dependency of the perturbed matrix elements on the value of #13, #23 and 
Qyi- We find that, although, (1, 1), (1, 2), and (1, 3) elements of the perturbed matrix do depend on 
all the three mixing angles, the #12 dependency is quite small compared to that of 6*13 and #23 . This 
is, in fact, true for normal as well as inverted hierarchy of the neutrino mass spectrum. However, the 
(2, 2), (2, 3), and (3, 3) elements depend only on the value of 623. It is expected as the perturbation 
to (2, 2), (2, 3), and (3, 3) elements of the tree level mass matrix was significantly smaller compared 
to (1, 1), (1, 2), and (1, 3) elements. Similar conclusions can be drawn for inverted hierarchy as well. 

We then use the results of bottom-up approach as guide-lines for determining models of neutrino 
masses and mixings. We introduce a typical mixing matrix Um = Utbm V for the neutrino mass 
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FIG. 8: Range of X(i,j) as a function of f3 for a = 0. The theory ranges are denoted by red whereas the 
experimental ranges from the bottom up approach are denoted by green. 

matrix, where V(a,f3, 7) is the perturbed mixing matrix. We fix the parameters a, (3, and 7 from 
the oscillation data obtained using the bottom up approach. We find that once we set a = 0, we 
do not get a good fit with the experimental data that are derived using the bottom up approach. 
Whereas, we can set (3 and 7 to be zero and still get a reasonably good fit. Since (3 and 7 can be 
set to zero, their ordering is not important. We further note that, if V is paramterized by more 
than one variable then the perturbations are arbitrary. In other words, those can not be treated as 
true perturbations with respect to their tree level mixing angles. 
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